Energy spectra of gamma-rays, electrons and neutrinos produced at interactions of 

relativistic protons with low energy radiation 



S.R. Kelnerf] and F.A. Aharoniaif] 

We derived simple analytical parametrizations for energy distributions of photons, electrons, and 
neutrinos produced in interactions of relativistic protons with an isotropic monochromatic radiation 
field. The results on photomeson processes are obtained using numerical simulations of proton- 
photon interactions based on the public available Monte-Carlo code SOPHIA. For calculations of 
energy spectra of electrons and positrons from the pair production (Bethe-Heitler) process we suggest 
a simple formalism based on the well-known differential cross-section of the process in the rest frame 
of the proton. The analytical presentations of energy distributions of photons and leptons provide 
a simple but accurate approach for calculations of broad-band energy spectra of gamma-rays and 
neutrinos in cosmic proton accelerators located in radiation dominated environments. 

PACS numbers: 12.20.Ds, 13.20.Cz, 13.60.-r, 13.85. Qk 



I. INTRODUCTION 

In astrophysical environments the density of low- 
energy radiation often exceeds the density of the gas 
component. Under such conditions, the interactions of 
ultrarelativistic protons and nuclei with radiation can 
dominate over interactions with the ambient gas. These 
interactions proceed through three channels: (i) inverse 
Compton scattering, jry —> jry' ', (ii) electron-positron pair 
production, p + 7 — > pe + e~ , and (hi) photomeson pro- 
duction, j»7 — * N + kn. While the inverse Compton scat- 
tering does not have a kinematic threshold, the electron- 
positron pair production and the photomeson produc- 
tion processes take place when the energy of the pho- 
ton in the rest frame of the projectile proton exceeds 
2m e c 2 ~ 1 MeV and m Tr c 2 (l + m 7r /2m p ) ~ 145MeV, re- 
spectively. 

The process of inverse Compton scattering of protons is 
identical to the inverse Compton scattering of electrons, 
but the energy loss-rate of protons is suppressed by a fac- 
tor of (m e /m p ) 4 w KT 13 . At energies above the thresh- 
old of production of electron-positron pairs this process 
is four orders of magnitude slower compare to the losses 
caused by pair-production. Therefore, generally the in- 
verse Compton scattering does not play significant role 
even in extremely dense radiation fields. 

The cross-section of (e + , e~ ) pair production (often 
called as Bethe-Heitler cross-section) is quite large, but 
only a small (< 2m e /m p ) fraction of the proton energy is 
converted to the secondary electrons. The cross-section 
of photomeson production is smaller, but instead a sub- 
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stantial (10 percent or more) fraction of the proton en- 
ergy is transferred to the secondary product. As a result, 
when the proton energy exceeds the 7r-meson production 
threshold, the hadronic interactions of protons dominate 
over the pair-production. 

The cross-section of pair-production is calculated with 
a very high accuracy using the standard routines of quan- 
tum electrodynamics. The cross-sections of photomeson 
processes are provided from accelerator experiments and 
phenomenological studies. Generally, for astrophysical 
applications the data obtained in fixed target experi- 
ments with gamma-ray beams of energies from 150 MeV 
to 10 GeV are sufficient, especially in the case of broad- 
band spectra of target photons, when the hadron-photon 
interactions are contributed mainly from the region not 
far from the energy threshold of the process. 

The energy losses of protons in the photon fields, in 
particular in the context of interactions of highest energy 
cosmic rays with 2.7 K cosmic microwave background ra- 
diation (CMBR), have been comprehensively studied by 
many authors (see e.g. Refs. (]], 0, p|, Less attention 

has been given to calculations of the products energy dis- 
tributions. Partly this can be explained by the fact that 
the cross-sections of secondary electrons and gamma-rays 
with the ambient photons significantly exceeds the cross- 
sections of interactions of protons with the same target 
photons. Therefore the electrons and gamma-rays cannot 
leave the active regions of their production, but rather 
trigger electromagnetic cascades in the surrounding ra- 
diation and magnetic fields. The spectra of gamma-rays 
formed during the cascade development are not sensitive 
to the initial energy distributions, and therefore simple 
approximate approaches (see e.g. Refs. ||) can provide 
adequate accuracies for calculations of the characteris- 
tics of optically thick sources. This does not concern, 
however, neutrinos which freely escape the source and 
thus have a undistorted imprint of parent protons. More- 
over, at some specific conditions, the secondary electrons 
from the pair- and photomeson production processes can 
cool mainly through synchrotron radiation for which the 
source can be optically thin (see e.g. J?], ||, ||, [h]]). Thus 
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the exact calculations of spectra of secondary electrons 
are quite important, since the synchrotron radiation of 
these electrons carry direct information about the energy 
spectra of accelerated protons. 

The interactions of hadrons with radiation fields can 
be effectively modeled by Monte Carlo simulations of 
characteristics of secondary products. In particular, the 
SOPHIA code provides an adequate tool for the com- 
prehensive study of high energy properties of hadronic 
sources in which interactions of ultrarelativistic protons 
with ambient low-energy photons dominate over other 
processes. At the same time, it is useful to have a com- 
plementary tool for study of radiation characteristics of 
hadronic sources, especially when one deals with sim- 
ple scenarios, e.g. interactions of protons with a homo- 
geneous and isotropic source of radiation in which the 
hadronic cascades (i.e. the next generation particles) do 
not play an important role. Motivated by this objective, 
in this paper we develop a simple approach based on the 
description of the energy distributions of final products 
from the photomeson and electron-positron pair produc- 
tion processes in analytical forms which can be easily in- 
tegrated in any complex model of hadronic interactions 
in high energy astrophysical sources. To a certain extent, 
this paper can be considered as a continuation of our first 
paper fll2|| where we obtained analytical presentations for 
proton-proton interactions. 

II. PHOTOMESON PROCESSES 

The formation of high energy gamma-rays, elec- 
trons and neutrinos in photomeson interactions proceeds 
through production and decay of non-stable secondary 
products, mainly 7r° and ^-mesons: 

7 + p — > n Q TT° + n + 7r + + n_TT~ + • • • (1) 

where n , n + and n_ are the numbers of produced pions. 
Hereafter we will assume that the density of the ambient 
medium is sufficiently low, thus the pions decay before 
they interact with the surrounding gas, radiation and 
magnetic fields. We will also assume that (i) both the 
relativistic protons and the target low-energy photons are 
isotropically distributed, and (ii) the energy of colliding 
particles 

e -C m^c 2 , E p > m p c 2 , (2) 

where e is the energy of the target photon, is the 
mass of 7r-meson (we will assume that mj = and 
E p and m p are the energy and the mass of proton. Those 
conditions, which are always satisfied in astronomical en- 
vironments, allows us to obtain simple analytical presen- 
tations for energy distributions of the final products of 
7r- decays - photons and leptons. 

The total cross-section of inelastic j-p interactions is 
a function of the scalar Sly = (k ■ p)/m p , where k and p 
are four-momenta of the photon and proton. The scalar 
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Figure 1: The total cross-section of the inelastic 7-p inter- 
actions as a function of energy of gamma-ray in the proton 
rest frame. The calculations have been performed using the 
routine of the code SOPHIA 0. 



S 1 is simply the energy of the photon in the proton rest 
frame. The total cross-section <t(S' 1 ), calculated using the 
routines of the code SOPHIA is presented in Fig. [I]. 
Although the production of pions dominate in the 7-p 
interactions, some other channels, in particular the ones 
leading to production of K- and 77-mesons, contribute 
noticeably (up to 10 to 20 %) to the overall production 
of photons and leptons. These channels are taken into 
account in our calculations presented below. They are 
based on the SOPHIA code which allows simulations of 
all important processes belongs to 7-p interactions. 



A. Production of gamma-rays 

The inclusive cross-section of production of 7r°-mesons 

da^o = S(k, p, p^o) d^p^o (3) 

depends on the momenta k, p, p„.o of the photon, proton 
and 7r°-meson, respectively. Let's denote by 

dw — W(p n o, p 7 ) d 3 p 7 (4) 

the probability of decay of a 7r°-meson with momentum 
p^o into a gamma-ray photon of momentum p 7 in the 
volume d 3 p 7 of the momentum space. Then 

da 1 = 2d 3 p 7 J S(k,p,p 7r o)W(p 7r o,p J ) d 3 p n o (5) 

can be treated as the inclusive cross-section of 7-rays 
production through the chain 7 + p — > n 7r° + ••■—» 
2n 7H . 

For isotropic distributions of initial particles, the final 
products will be isotropically distributed as well. There- 
fore for the determination of the energy spectra of 7-rays 
we can use the inclusive cross-section given by Eq. (ph but 
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integrated over the gamma-rays emission angles: 

dE 

da-y = Gry(E p , e, cos 9, E y ) —- 1 

E n 



(6) 



The function G 7 depends on E p , e, E 1 , and the angle 
9 between the momenta of colliding proton and photon. 
The corresponding differential interaction rate is 



dw. 



dE^ 



~ — c (1 — cos 9) Gj(E p , e, cos 9, E 7 ) — - 

E„ 



(7) 



Let fp(Ep) and / P h(e) be functions characterizing the 
energy distributions of initial protons and photons, i.e. 
fp(Ep) dE p and / P h(e) de are the numbers of protons and 
photons per 1 cm 3 in the energy intervals dE p and de. 
respectively. Since it is assumed that the target photons 
are isotropically distributed, their angular distribution is 
described as dVt/Ait. Then the production rate of 7-rays 
(i.e. number of 7-rays in the energy interval {E~,E 1 + 
dE-y) per sec, per 1 cm 3 ) can be obtained after integration 
of Eq. (0) over energies of protons and target photons, as 
well as over the solid angle effi: 



dN 7 (E 7 ) = dE,J ^^de fp{E p ) / ph ( 

c(l - cos6) G y (E p , e, cos9 7 E 7 ) . (8) 
Let's introduce the function $ 7 defined as 

c (1 - cos#) G~{E V , e,cos#,i? 7 ) — , (9) 



where 



4eE p 
m 2 c 4 



E~. t 
E„ 



Then Eq.(ph can be written in the following form 



dN 7 

HeZ 



/dE 
fp(E p ) /ph(e) $ 7 (ry, x)-^de. 



(10) 



(11) 



Note that $ 7 can be treated as a function of two (but not 
three) variables. This is connected with the fact that, as 
we assume, at E p — > 00 and e — + 0, E 1 — * 00, and for 
fixed x and 77, the function <i> 7 should have a certain limit. 
In other words, at the rescaling 



Ep — > XEp Ej — > XEj , e — > e/A , 



(12) 



where A is an arbitrary number (A > 1), the function 
<I> 7 is not changed, therefore the following relation takes 
place. 

/(l-cos0) G 7 (A£ p ,e/A,cos6>,A£ 7 )^ = 
= f(l-cos9)G^(Ep,e,cos9,E^)dQ. (13) 

Numerical calculations based on the code SOPHIA |ll[] 
show that already at 

e < 1(T 3 m^c 2 , Ep > 10 3 m p c 2 (14) 



the relation given by Eq.|TJ) is readily fulfilled. 

Some features of the energy spectra of 7-rays can be 
understood from the analysis of the kinematics of the pro- 
cess. In particular, for production of a single 7r°-meson 
the following condition should be satisfied 



2 eE p (1 — f3 p cos 0) > (2 to w to p + m 2 ) c 4 , 



(15) 



where (3 p - is the proton speed (in the units of c). 

Since the protons are ultrarelativistic, we will assume 
(3 P = 1. If the condition of Eq.( [15|) is not satisfied, 
then the interaction rate given by Eq. (J7|) is equal to zero. 
Therefore, for the case of 4,eE p < (2m 7r m p + m 2 )c 4 , the 
function <I> 7 = 0. The integration of Eq. (pi]) should be 
performed over the region 1 



V > Vo = 2 • 



0.313. 



(16) 



As it follows from kinematics of production of a single 
pion, the energy of the latter appears within 



-^7T mill ^ -^7T — Eft max , 



where 



E-Kl 



— Er, Xj 



En min — E p X _ 



(17) 



(18) 



are the maximum and minimum energies respectively, 
with 

1 



2(1 + 77) 



rj + r 2 ± \/ (r\ — r 2 — 2r) (rj — r 2 + 2r) , 

(19) 

where r = m n /m p m 0.146. 

Let's consider now the general case when the total 
mass of particles produced in proton-photon interactions 
is M + m^. For example, the single-pion production im- 
plies M — m p , while in the case of two-pion production 
M = trip + 771^, etc. In this case the maximum and min- 
imum energies of pions are 



EL 



with 



Er, x_i 



1 

2(1 + »?) 



EL 



Ep x_ 



i-ir ± 



(20) 



± V(J7 + 1 - (R + r) 2 ) (n + l-(R- r) 2 ) 



(21) 



where R — M/m p . In particular, for R = 1 Eqs. jl9| ) 
and (^lj) coincide. For R > 1, we have the following 
inequalities 



Ejr 



Ejr 



(22) 



1 For determination of the region allowed by kinematics, we will 
assume m 7T ± = m n o = 0.137GeV. Since on the boarder of this 
region the function "I^ = 0, this approximation does not affect 
the accuracy of numerical calculations. 
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The decay of ultrarelativistic 7r°-mesons with energy 
distribution J„(E„) within the limits given by Eq.(17), 
results in the energy spectrum of 7-rays 



dN 1 
dE^ 

where 

Ei = max E-y , ■ 



= 2 



Ejr 



J n (E n ) 



(23) 



£1 



! E n I 



EL < En 1 



(24) 



Below we will be interested in 7-rays with energy E 1 > 
m-n/2. While in the energy r^ngc E*y <C E^ min 

the diffe- 
rential spectrum of 7-rays is flat, 



dN 1 
dEL 



"r*dE„ 
2 / -g— J^Eir) 



(25) 



within the interval E nm \ n < E-y < -E^max the spectrum 
decreases with energy, 



d/V 7 
dEL 



= 2 



dE n 

Ejr 



Jn(E n ) 



(26) 



The above quantitative features of the spectrum of 7- 
rays appear quite useful for the choice of approximate 
analytical presentations. The results of numerical cal- 
culations of the function $ 7 based on simulations using 
the code SOPHIA can be approximated, with an ac- 
curacy better than 10% by simple analytical formulae. 
Namely, in the range x~ < x < x + 



<& 1 {rj 1 x) = i? 7 exp<j 









| S 7 


Hi:)'. 









In 



i + ir 



2.5+0.4 In(7j/r) ) 



where 



(27) 



(28) 



At low energies, x < sc_, the spectrum does not depend 
on x, 



$ 7 (77,21) = B 1 (in 2) 



2.5+0.4 ln(r;/%) 



(29) 



Finally in the range x > x+ the function $ 7 = 0. 

All three parameters i? 7 , s 7 and <5 7 used in this presen- 
tation are functions of 77. The numerical values of these 
parameters are shown in Table [i]. At 77/770 — 1, i.e. at 
the threshold of 7r°-meson production, B 1 = 0. In Fig. || 
we show the functions x^-y(rj,x) obtained with the code 



Table I: Numerical values of parameters 73-, , s 7 and <5 7 char- 
acterizing the 7-ray spectra given by Eq.dTtl). 



V/Vo 


s 7 


S~, 


B 1 , 


cm 3 


h 


1.1 


0.0768 


0.544 


2.86 


• 10" 


-19 


1.2 


0.106 


0.540 


2.24 


■ 10- 


-18 


1.3 


0.182 


0.750 


5.61 


• 10" 


-18 


1.4 


0.201 


0.791 


1.02 


• 10- 


-17 


1.5 


0.219 


0.788 


1.60 


■ 10" 


-17 


1.6 


0.216 


0.831 


2.23 


• 10" 


-17 


1.7 


0.233 


0.839 


3.10 


• 10" 


-17 


1.8 


0.233 


0.825 


4.07 


• 10" 


17 


1.9 


0.248 


0.805 


5.30 


• 10" 


-17 


2.0 


0.244 


0.779 


6.74 


• 10" 


-17 


3.0 


0.188 


1.23 


1.51 


• 10" 


-16 


4.0 


0.131 


1.82 


1.24 


• 10" 


-16 


5.0 


0.120 


2.05 


1.37 


• 10" 


-16 


6.0 


0.107 


2.19 


1.62 


■ 10" 


-16 


7.0 


0.102 


2.23 


1.71 


• 10" 


-16 


8.0 


0.0932 


2.29 


1.78 


• 10" 


-16 


9.0 


0.0838 


2.37 


1.84 


• 10" 


-16 


10.0 


0.0761 


2.43 


1.93 


• 10" 


-16 


20.0 


0.107 


2.27 


4.74 


■ 10" 


-16 


30.0 


0.0928 


2.33 


7.70 


• 10" 


-16 


40.0 


0.0772 


2.42 


1.06 


• 10" 


-15 


100.0 


0.0479 


2.59 


2.73 


• 10" 


15 



SOPHIA (histograms) and the analytical presentations 
given by Eqs. ( p7j ) and (^) for two values of 77. 

Eq.([ll]) provides a simple approach for calculations of 
7-ray spectra for arbitrary energy distributions of ultra- 
relativistic protons and ambient photons. The param- 
eters -B 7 , s 7 and S-y are quite smooth functions of 77, 
thus for calculations of these parameters at intermedi- 
ate values of 77 one can use linear interpolations of the 
numerical results presented in Table |j| Note that for in- 
teractions of protons with 2.7 K CMBR, the results pre- 
sented in Table | allow calculations of 7-ray spectra up 
to - 10 21 eV. 



B. Production of electrons and neutrinos 

The production of leptons in proton-photon interac- 
tions is dominated by the decay of secondary charged 
pions. In analogy with Eq. (p"T|) , the spectrum of each 
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X *y(.V.X) 



X *y(.V.X) 





Figure 2: Gamma-ray spectra produced in photomeson interactions calculated for two values of r\ = 4eE p /rripC 4 . Solid lines are 
calculated using the analytical presentations given by Eq.rtlll), the histograms are from simulations using the SOPHIA code. 



type of leptons can be presented in the form 
N l {E l )dE l = 
= dE l J f p (E p ) / ph (e) $ ; (r,,x)^de, (30) 

where 77 is determined in Eq.(|l(i|), x = E l jE v and I im- 
plies one of the following symbols: e + , e~, v el v e . 

As in the case of 7-rays, the energy range of leptons, 
for the fixed values of E p and e is determined by kinemat- 
ics. We obtained analytical forms for the energy spectra 
for all lepton types, using numerical results of simu- 
lations the SOPHIA code. In the range of x'_ < x < x' + 



$i(r],x) = B l exp< -s l 



In 



where 



x — x 



(31) 



(32) 



In range of x < x'_ , the function does not depend on 
x: 



<S> l { V ,x) = B l (In 2) 



(33) 



For x > x' + the function = 0. The analytical presen- 
tation in the form of Eq.(|3l|) contains three parameters 
s,, 5r and B l which themselves are functions of rj. The 
numerical values of these parameters obtained with the 
method of least squares are tabulated in Table O (for e + , 



1. Energy spectra of positrons, muon antineutrinos and 
electrons neutrinos 



For e + , 
form 



with 



v and v e , the parameter ip is presented in the 



ip = 2.5 + 1.4 In 



- 



x'_ = and x' + = x + , 



(34) 
(35) 



where x+ and x_ are determined from Eq.([l9|); r/ is 
defined in Eq. (|i"6|) . Note that here tp is different than 
the relevant function in Eq.(||). In Fig. | and Fig. | the 
analytical presentations of distributions x $ e + and x $ p 
are compared with Monte-Carlo simulations based on the 
SOPHIA code. 

In the range 77 < 4 ^- + 4 fa^ = 2.14 770, only a sin- 
gle 7r + -meson can be produced. It decays to tt + — > [i^v^. 
The positrons and muon antineutrinos are produced from 
the decay /_t + — > e + u^v e . Since the spectra of e + and 
from the decay of coincide (see e.g. [|l3]), the param- 
eters in the Tables || calculated for e + and for small 
77 should be identical. The slight difference at 77 < 2?yo 
is explained by fluctuations related to the the statistical 
character of simulations. 

At 77 > 2.14t7o, a new channel is opened for produc- 
tion of v^, because of production of 7r~-mesons and their 
decay: it~ —> ^TV^. Therefore for large values of 77 the 
parameters characterizing e + and differ significantly 
(see Table ||. This can be seen from comparison of re- 
sults presented in right panels of Fig. || and |; for r\ = 30, 
x $ p l significantly exceeds x <& e + . 



2. Muon neutrinos 



and v e ) and Table ||| (for e~ and v e ). The values The distribution for is described by Eq.Q) with the 



of a;± and tp are given below for each type of leptons. 



same function ip as for e + , v and v e , given by Eq.(p4| 
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Table II: Numerical values of parameters s h 5 t , and B l for e + , v^, v e . 



V/Vo 


s e + 




B e+ , 


cm 3 /s 






Bd^ , 


cm 3 /s 








cm 3 /s 






B V<1 , 


cm 




1.1 


0.367 


3.12 


8.09 


■ 10" 


19 


0.365 


3.09 


8.09 


■ 10" 


19 


0.0 


0.0 


1.08 


■ 10" 


18 


0.768 


2.49 


9.43 


10 


-19 


1.2 


0.282 


2.96 


7.70 


■ 10" 


18 


0.287 


2.96 


7.70 


• 10" 


18 


0.0778 0.306 


9.91 


■ 10" 


18 


0.569 


2.35 


9.22 


10 


-18 


1.3 


0.260 


2.83 


2.05 


■ 10" 


17 


0.250 


2.89 


1.99 


■ 10" 


17 


0.242 


0.792 


2.47 


• 10" 


17 


0.491 


2.41 


2.35 


10 


-17 


1.4 


0.239 


2.76 


3.66 


■ 10" 


17 


0.238 


2.76 


3.62 


• 10" 


17 


0.377 


1.09 


4.43 


■ 10" 


17 


0.395 


2.45 


4.20 


10 


-17 


1.5 


0.224 


2.69 


5.48 


■ 10" 


17 


0.220 


2.71 


5.39 


• 10" 


17 


0.440 


1.06 


6.70 


■ 10" 


17 


0.31 


2.45 


6.26 


10 


-17 


1.6 


0.207 


2.66 


7.39 


■ 10" 


17 


0.206 


2.67 


7.39 


■ 10" 


17 


0.450 


0.953 


9.04 


• 10" 


17 


0.323 


2.43 


8.57 


10 


-17 


1.7 


0.198 


2.62 


9.52 


■ 10" 


17 


0.197 


2.62 


9.48 


• 10" 


17 


0.461 


0.956 


1.18 


• 10" 


16 


0.305 


2.40 


1.13 


10 


-16 


1.8 


0.193 


2.56 


1.20 


■ 10" 


16 


0.193 


2.56 


1.20 


■ 10" 


16 


0.451 


0.922 


1.32 


• 10" 


16 


0.285 


2.39 


1.39 


10 


-16 


1.9 


0.187 


2.52 


1.47 


■ 10" 


16 


0.187 


2.52 


1.47 


• 10" 


16 


0.464 


0.912 


1.77 


■ 10" 


16 


0.270 


2.37 


1.70 


10 


-16 


2.0 


0.181 


2.49 


1.75 


■ 10" 


16 


0.178 


2.51 


1.74 


■ 10" 


16 


0.446 


0.940 


2.11 


• 10" 


16 


0.259 


2.35 


2.05 


10 


-16 


3.0 


0.122 


2.48 


3.31 


■ 10" 


16 


0.123 


2.48 


3.38 


■ 10" 


16 


0.366 


1.49 


3.83 


■ 10" 


16 


0.158 


2.42 


3.81 


10 


-16 


4.0 


0.106 


2.50 


4.16 


■ 10" 


16 


0.106 


2.56 


5.17 


■ 10" 


16 


0.249 


2.03 


5.09 


■ 10" 


16 


0.129 


2.46 


4.74 


10 


-16 


5.0 


0.0983 


2.46 


5.57 


■ 10" 


16 


0.0944 


2.57 


7.61 


■ 10" 


16 


0.204 


2.18 


7.26 


■ 10" 


16 


0.113 


2.45 


6.30 


10 


-16 


6.0 


0.0875 


2.46 


6.78 


■ 10" 


16 


0.0829 


2.58 


9.57 


• 10" 


16 


0.174 


2.24 


9.26 


■ 10" 


16 


0.0996 


2.46 


7.65 


10 


-16 


7.0 


0.0830 


2.44 


7.65 


• 10" 


16 


0.0801 


2.54 


1.11 


■ 10" 


15 


0.156 


2.28 


1.07 


■ 10" 


15 


0.0921 


2.46 


8.61 


10 


-16 


8.0 


0.0783 


2.44 


8.52 


■ 10" 


16 


0.0752 


2.53 


1.25 


• 10" 


15 


0.140 


2.32 


1.19 


■ 10" 


15 


0.0861 


2.45 


9.61 


10 


-16 


9.0 


0.0735 


2.45 


9.17 


■ 10" 


16 


0.0680 


2.56 


1.36 


■ 10" 


15 


0.121 


2.39 


1.29 


■ 10" 


15 


0.0800 


2.47 


1.03 


10 


-15 


10.0 


0.0644 


2.50 


9.57 


■ 10" 


16 


0.0615 


2.60 


1.46 


■ 10" 


15 


0.107 


2.46 


1.40 


■ 10" 


15 


0.0723 


2.51 


1.10 


10 


-15 


30.0 


0.0333 


2.77 


3.07 


■ 10" 


15 


0.0361 


2.78 


5.87 


• 10" 


15 


0.0705 


2.53 


5.65 


■ 10" 


15 


0.0411 


2.70 


3.55 


10 


-15 


100.0 


0.0224 


2.86 


1.58 


• 10" 


14 


0.0228 


2.88 


3.10 


• 10" 


14 


0.0463 


2.62 


3.01 


■ 10" 


14 


0.0283 


2.77 


1.86 


10 


-14 



x * e +(7),a;) x ^ e *{-r),x) 




Figure 3: Energy spectra of positrons produced in photomeson interactions calculated for two values of rj — AeEp/rripC 4 . Solid 
lines are calculated using the analytical presentations, the histograms are from Monte-Carlo simulations using the SOPHIA 
code. 



but with different parameters x'±: 

' 0.427 x+, p<2.U, 
(0.427 + 0.0729 (p - 2.14)) x+ , 2.14 < p < 10 , 
x+ , p > 10 , 

(36) 



x+ = < 



where p = 77/770, and 



0.427 x_ 



(37) 



The difference of values of x± for and v ^ appears for 
the following reason. At the decay ir + — ► p + v^ the maxi- 
mum energy of £> M is equal to (1— m^/m^.) w 0.427^, 
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where is the mass of muon. On the other hand, the 
maximum energy of v^, produced at decay of the muon 
is equal to E w . With an increase of the parameter 77, tt~- 
mesons start to be produced, the decay of which leads 
to with maximum energy comparable to the energy of 
the pion (a detailed discussion of these questions can be 
found in |12j). 

In Fig. g the function x $„ (77, x) is shown. The 
histograms are from Monte-Carlo simulations using the 
SOPHIA code, and the solid lines correspond to the an- 
alytical approximations. 



3. Electrons and electron antineutrinos 

The electrons and electron antineutrinos are produced 
through the decay p~ , which in its turn is a product 
of the decay of 7r _ -meson. Therefore, for production of 
e~ and P e , at least two pions should be produced. The 
production of two pions is energetically allowed if 

»7>< in = 4r(l + r)w2.14r to . (38) 



The maximum and minimum energies of the pion corre- 
spondingly are 

E-k max — ^max : -^tt min — *^min ^pi > (*^) 

where 

= 2(1 + 77) ( ?/ ~ 2 r ± V^FMT+7))) . (40) 

Then x' + = x max and x'_ = x m j n /2. Eq.([i"o|) is obtained 
from Eq.(pTl) if one sets R = 1 + r. These functions 
together with 

^ = 6fl-e 1 ' 5(4 - p)N ) 6(p-4), p=—, (41) 

determine the distributions for e~ and v e given in a gen- 
eral form by Eqs. (^Tj) - ([33]); O(p) is the Heaviside func- 
tion (B(p) = if p < and O(p) = 1 if p > 0). 
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Table III: Numerical values of parameters s t , S l and B t for 
electrons and electron antineutrinos 



V/Vo 


s e - 


8 e - 


B e - , 


cm 3 










cm 3 




3.0 


0.658 


3.09 


6.43 


• io- 


-19 


0.985 


2.63 


6.61 


io- 


19 


4.0 


0.348 


2.81 


9.91 


• io- 


-18 


0.378 


2.98 


9.74 


io- 


18 


5.0 


0.286 


2.39 


1.24 


■ 10" 


-16 


0.31 


2.31 


1.34 


io- 


16 


6.0 


0.256 


2.27 


2.67 


■ 10" 


-16 


0.327 


2.11 


2.91 


10" 


16 


7.0 


0.258 


2.13 


3.50 


■ 10" 


-16 


0.308 


2.03 


3.81 


io- 


16 


8.0 


0.220 


2.20 


4.03 


■ 10" 


-16 


0.292 


1.98 


4.48 


io- 


16 


9.0 


0.217 


2.13 


4.48 


• 10" 


-16 


0.260 


2.02 


4.83 


io- 


16 


10.0 


0.192 


2.19 


4.78 


• 10" 


-16 


0.233 


2.07 


5.13 


io- 


16 


30.0 


0.125 


2.27 


1.64 


• 10" 


-15 


0.135 


2.24 


1.75 


10" 


15 


100.0 


0.0507 


2.63 


4.52 


■ 10" 


-15 


0.0770 


2.40 


5.48 


io- 


15 



III. PHOTONS AND LEPTONS PRODUCED AT 
INTERACTIONS OF PROTONS WITH 2.7K 
CMBR 

In this section we compare the energy spectra of 
gamma-rays, neutrinos, and electrons produced at pho- 
tomeson interactions. For monoenergetic protons inter- 
acting with a radiation field with energy distribution 
/ph(e), the energy spectra of photons and leptons can be 
reduced to the calculation of a one-dimensional integral 



oo 

= f Zph(e) $(»7, x) de. 



(42) 



where e = r] m p c r/(4B ), and x = E/E p is the fraction 
of energy of the protons transferred to the given type of 
secondary particle; $ is one of the functions described in 
the previous section. 

In Fig. ^| and ^ we show the energy spectra of 7- 
rays and electrons (left panels) and all neutrino types 
(right panels) produced by protons of energy 10 20 eV and 
10 21 eV interacting with blackbody radiation of temper- 
ature T = 2.7 K. The results depend only on the prod- 
uct E p x T, therefore they can be easily rescaled to a 
blackbody radiation of an arbitrary temperature. The 
chosen radiation field and proton energies are of great 
practical interest in the context of origin and intergalac- 
tic propagation of ultrahigh energy cosmic rays. Because 
of interactions with the intergalactic radiation fields, ul- 
trahigh energy 7-rays achieve the observer from distances 
less than 1 Mpc (see e.g. Ref. 14 ). The electrons rapidly 



cool via synchrotron radiation or, in the case of very small 
intergalactic magnetic field, initiate electromagnetic cas- 
cades supported by interactions of electrons and gamma- 
rays with the 2.7 K CMBR. Only neutrinos freely pene- 
trate through intergalactic radiation and magnetic fields 
and thus carry a clear imprint of parent protons. 

In Fig. |8| we show the average number of secondaries 



(multiplicity) produced in one inelastic interaction of pro- 
tons with 2.7 K CMBR as a function of proton energy. 
The results of numerical calculations are obtained using 
the energy spectra of secondary photons, electrons and 
neutrinos, and the total cross-section shown in Fig ^. 
Note that below the threshold of production of two pions 
one should have the following relations 

- (n 7 ) + (n e +) = 1, (n e +) = (n % ) = (n Uft ) . (43) 

The results of calculations based on approximate analy- 
tical presentations of functions <I>; satisfy these relations 
with an accuracy of better than 5%. 

Note that at very low energies the average number 
of 7-rays (n 7 ) appears smaller than average number of 
positrons (n e +). This, at first glance unexpected result 
is actually a direct consequence of the experimental fact 
that near the threshold the total cross-section a 7r + of ir + 
production significantly exceeds a^o (see Fig. ||). 



IV. PRODUCTION OF ELECTRON-POSITRON 
PAIRS 

At energies below the photomeson production, the 
main channel of inelastic interactions for protons with 
ambient photons proceeds through the direct production 
of electron-positron pairs. In the rest frame of the proton, 
this process is described by the so-called Bethe-Heitler 
cross-section. In astrophysical environments, the process 
is more often realized when ultrarelativistic protons col- 
lides with low-energy photons, 



p + 7^e + + e +p 
The process is energetically allowed when 

7 p e > m e c 2 , 



(44) 



(45) 



where j p = E p /m p c 2 is the proton Lorentz factor, e is the 
soft photon energy, and m e is the mass of electron. The 
maximum energy of the electron (positron) is determined 
by the kinematics of the process 



1p 



1 + 47 p e/(m p c 2 ) 



1p ( 



(46) 

This equation is valid when 7 p 3> 1 and e ^ m p j p c 2 . In 
the interval 



m e c 2 7 p e <^ m p c 2 



the maximum electron energy is 



E„ 



This applies for E emax In the limit of 7 p e 3> 

Eemax — TYlpC 7p — Ep , 



(47) 



(48) 



(49) 
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x dN/dx , s 




x=E/E p 



Figure 7: The same as in Fig. ^[ but 



x dN/dx , s 




x=E/E p 



for a proton of energy E p = 10 eV 



i.e. the whole energy of the proton is transferred to one 
of the electrons. 

Let's denote by da the differential cross-section of the 
process. The interaction rate is 



dw 



„3 ( k ' P) , 



{k ■ Up 
Hp 



da , 



(50) 



where k and p are four-momenta of the photon and 
proton, Up — p/uipC is the four-velocity of the proton, 
(k-p) — eEp/c 2 — kp is the scalar product of four- vectors. 
Let assume that in a unit volume we have / p h(e) de dtt/Air 
photons between the energy interval (e, e -I- de) and mov- 
ing within the solid angle dil. Then the number of inter- 
actions per unit of time is 



V = ^/de^/ ph (e) 



(k ■ u p ) 
£7p 



da . 



(51) 



where the integration is performed over all variables. 

Below we perform calculations based on the following 
approach. If we are interested in a distribution of some 
variable £, which is a function ip of particle momenta, 



this distribution can be found introducing an additional 
8- function under the integral in Eq.(J5T|): 



dN 2 r dn (k ■ up) 

— = c de — f ph (e) ^ 

a£ J 47r £7 P 



8(£ -cp)da. (52) 



In particular, the energy distribution of electrons in the 
laboratory frame can be calculated using the following 
formula 



dN 
dE~, 



dn (k ■ Up) 
de — /ph(e) — 

47T 67 p 



8(E e -c(u lf -p e ))da, 
(53) 

where u« is the four-velocity of the laboratory frame, 
and p e is four-momentum of electron since the scalar 
c(m ; j -p e ) is equal to the energy of electron in the labora- 
tory system. The proton Lorentz-factor in the laboratory 
system also can be considered as a relativistic invariant: 
7 P = (u v ■ Up). 



if 

Note that the integral 



5 = / 8(E e - c(u lf ■ p e ))da 



(54) 



10 



6- , the result can be written in the form 



3 
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Figure 8: The multiplicity of photons and leptons produced 
in one interaction of a relativistic proton with 2.7 CMBR 



0.3 



0.2 



0.1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 

S y , GeV 

Figure 9: The total cross-sections of production of ir and 
7r°-mesons as a function of energy of incident gamma-ray in 
the rest frame of proton. The experimental points are taken 
from the site http://wwwppds.ihep. su:8001. 



is a relativistic invariant, so it can be calculated in any 
frame of coordinates. The differential cross-section da 
can be written in the simplest form in the rest-frame of 
the proton, therefore for calculations of S we will use this 
system of coordinate where 

c(u lf -p e ) = 7 P (£L - V p p- cos#_) . (55) 

Here E_ is the energy and p_ = »J E 2 _j c 2 — m 2 c 2 is the 

momentum modulus of electron in the rest frame of the 
proton, is the angle between the momenta of photon 
and electron. Therefore 

S = J da5(E e - 7p (£_ - V p p- cos<9_)) . (56) 
After integration over all variables, except for £L and 



S 



5{E e - lp {E_ 



W(w, J B_,cos6»_)x 

V p p- cos 9- ) ) dE_ d(cos 0_ ) 



where 



W(u),E_,cos6^) = 



d 2 a 



dE_ d(cos6L) 



(57) 



(58) 



is the double-differential cross-section as a function of en- 
ergy and emission angle of the electron in the rest frame 
of the proton; u> = (u p ■ k)/(m e c) is the energy of the 
photon the rest frame of proton in units m e c 2 . The func- 
tion W has been derived in the Born approximation in 
Refs. [[l5[ |l(|. The approach used in these papers de- 
scribes the production of an electron-positron pair by a 
photon in the Coulomb field which formally corresponds 
to the limit m p — > oo. However, we warn the reader 
that there is a misprint in the cross-section published in 
these papers, therefore we advise to use Eq.(lO) of paper 
by Blumenthal M , where the typo is fixed. Note that 
in Refs. g, |l5|, the system of units is used in which 
c = h = m e = 1. Since here we cite to certain equa- 
tions of these papers, in this section, in order to avoid a 
confusion, we use the same system of units. 

The presence of the 5-function in the integrand allows 
the integration over the variable d(cos0_), which gives 



S = 



1 



IpVp 



dE_ 



W(cj,E_,C), 



where 



1 P E- 



E„ 



IpVpP- 



(59) 



(60) 



After substituting Eq.(|59|) into Eq. (|53|) , and using the 
relation 



uj = (u p ■ k) = ej p (1 — cos ( 



(61) 



it is convenient to perform the integration over to instead 
of integration over the angle. Then, for ultrarelativistic 
protons (7p 3> 1), we obtain 



dN _ 1 
dEe~2ll 



dc 



/ph(e) 



27 P £ 



dujux 



(•r „+ger 



dE_ 



(62) 



When substituting Eq. (pOj) into ( |62| ) we set V p = 1, and 
correspondingly £ = (7 P -EL — i? e )/(7 p p_). The integra- 
tion limits in Eq. ( |62| ) are found from the analysis of kine- 
matics. 
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In the case of monoenergetic target photon field, 

f ph (e') = CS(e'-e) (63) 

the energy distribution of electrons can be written in the 
form of the double-integral 



E 2 dN/dE , eV s" 1 



dN _ C 
dEe ~ 2^p 



27 p e 



du> u> x 



dE„ 
P- 



with the following kinematic condition 
tefiE e > (% + E e ) 2 . 



(64) 



For the important case of Planckian distribution of tar- 
get photons, 



/ph(e) 



1 



7T 2 e £ / feT - 1 ' 



(65) 



the expression can be simplified. Indeed, rewriting the 
first term in the integrand of Eq. ( |62| ) in the form 



i(l - e- £ / fcT ) 



(66) 



we can perform integration over de by parts, which after 
simple transformation leads to the energy spectrum of 
electrons in the form of double-integral 



dN 
dK 



kT 



duj lu In ( 1 — e 



-w/(2 7p fcT) x 



dE_ 



W(u,E_,0. 



(67) 



In the Born approximation used in g, [t(|, the en- 
ergy and angular distributions of electrons and positrons 
are identical, therefore Eqs. ( |62"|) and (|67]) do not distin- 
guish between electrons and positrons. Let's discus the 
condition of applicability of Eqs. ( |62| ) and (|67|). In the 
proton rest system, the cross-section of production of an 
electron-positron pair by the proton and in the Coulomb 
potential coincide for all emission angles of pairs, when 
u> <C uip. In the laboratory frame this is equivalent to the 



condition e"f p <C m p , which, taking into account Eq.(48) 



can be written in the form E enls , x <C E p . Thus, the above 
obtained results can be applied to production of electrons 
and positrons when E e E p . 

In Fig. [l^ we show the energy distributions of electrons 
and positrons produced in interactions of protons of three 
different energies with the 2.7 K CMBR: 6.4 x 10 19 eV, 
10 20 eV and 3 x 10 20 eV. Note that the energy of pri- 
mary proton 6.4 x 10 19 eV is interesting in the sense 
that at this energy the loss-rates of protons, E~ x dE p j 'dt, 



10' 



10" 




10' 



Figure 10: Energy distributions of electrons and positrons 
(iV e = N+ +N-) produced at interactions of protons with 2.7 
K CMBR. Dashed lines and dot-dashed lines correspond to 
the pair production and photomeson processes, respectively. 
Solid lines show the sums of two contributions. The curves 
correspond to three energies of protons: 6.4 x 10 19 eV, 10 20 eV 
and 3 x 10 20 eV. 



due to pair production and photomeson production are 
equal. The spectral energy distribution of electron and 
positrons from the pair production process, E 2 dN/dE, 
has a bell-type shape with a broad maximum around 
(m e /m p )R p ~ 10~ 3 E P . This spectrum is quite different 
from the dN/dE oc E' 1 ^ type energy dependence as it 
was hypothesized in |l(J . Fig. [it] demonstrates that while 
the low energy range of electrons (positrons) is dominated 
by the process of pair production, at higher energies 
the main contribution comes from photomeson processes. 
Fortunately, in the energy range where Eqs. ([32]) and ( |67| ) 
are not valid, the contribution of pair-production to the 
spectrum of electrons is negligible compared to the con- 
tribution of photomeson processes. 

It should be noted that the ^-functional approxima- 
tion, which is often used for qualitative estimates of char- 
acteristics of products of high energy interactions, in this 
specific process does not provide adequate accuracy. The 
reason is that the electrons produced in a single act of in- 
teraction have very broad energy distribution. The calcu- 
lations show that the ^-functional approximation leads to 
significant deviation from exact result given by Eq. (|62[) , 
even when one takes into account the energy dependence 
of the average fraction of the proton energy transferred 
to the electron. 



A. Energy losses 

The analytical presentations of the energy spectra of 
stable products of interactions of protons with ambient 
low energy photons allows us to calculate the energy 
losses of protons in a radiation field with arbitrary en- 
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\dE p /dt\/E p , yr" 1 




E p . eV 



Figure 11: The average energy lose rates of protons in the 
CMBR with temperature 2.726 K. The lines are obtained us- 
ing our method of integration of energy spectra of all final 
(stable) secondaries, the points are from H. They are ob- 
tained from Monte Carlo simulations of interactions of pro- 
tons with the CMBR photons using the SOPHIA code. The 
dashed and dash-dotted lines describe the energy losses due 
to pair-production and photomeson interactions, respectively, 
the solid lines represent the sum of these two contributions. 



dN/dt , yr -1 (x) 




Figure 12: The interaction rates of protons with photons 
of 2.7 K CMBR (left axis) and the coefficient of inelasticity 
(right axis). The curves la and 2a are the electron-positron 
and photomeson production rates, respectively. The curves 
lb and 2b are the average energy lost by a proton of given 
energy due to pair production and photomeson production, 
respectively. 



ergy distribution, 



1 



dEr, 



dt 



1 CO 

J dx x J de / p h(e) $(r/, x) 



(68) 



where e m in = 



/(4Ep), and ^(rj, x) is the sum of all 



seven energy distribution (relevant to 7, e + , e~, v^, v^, 
u e , i7 e ) derived in Sec. 2. Eq.(^) describes the average en- 
ergy losses transferred to gamma-rays and leptons. In or- 
der to calculate the energy losses due to pair-production 
one should multiply Eq.(67) to 2E e and integrate over 
dE e . 

Calculations of energy losses of protons can be per- 
formed directly, without intermediate calculations of en- 
ergy distributions of secondary photons and leptons. In 
this regard, the energy losses of protons in 2.7 K CMBR 
have been studied in great details by many authors, in 
particular by Berezinsky and coauthors based on a 
semi-analytical method of calculations and Stanev et al. 
H based on Monte Carlo simulations using the SOPHIA 
code. Therefore it is interesting to compare our results 
with direct calculations of energy losses. In Fig. [n] we 
show the energy lose rate of protons in the blackbody 
radiation field with temperature T = 2.726 K. For com- 
parison, we show the result of calculations performed us- 
ing the code SOPHIA Q. The agreement of two calcu- 
lations is an indirect test of a good accuracy of the ob- 
tained above approximate analytical presentations for en- 
ergy distributions of stable products from proton-photon 
interaction. 



In Fig. [jj] we show the interaction rate of protons 
with 2.7 K CMBR, as well as the fraction of energy lost 
by the proton per interaction (the so-called inelasticity 
coefficient). Close to the threshold of pair production 
around E ~ 10 18 eV, (x) e+e - = 2m e /m p « 1.1 • 1(T 3 , 
as it expected from the kinematics of the process. How- 
ever, with an increase of energy, (x) v gradually decreases 
down to 1(T 4 at 10 20 eV. This effect has been noticed 
also in [lTj. In the case of photomeson production 
the inelasticity coefficient has a quite different behav- 
ior. At the threshold, (a;)^ increases from the value of 
m„/(m p + m n ) « 0.13 to approximately 0.4 at energy 
10 22 eV. Therefore, despite the fact that the cross-section 
of pair production significantly, by two orders of magni- 
tude, exceeds the cross-section of photomeson produc- 
tion, the energy losses at high energies are dominated by 
photomeson interactions. 



CALCULATIONS FOR POWER-LAW 
DISTRIBUTION OF PROTONS 



Instead of integrating Eqs. ( |ll|) and (|30|) over de, it 
is more convenient to perform integration of these equa- 
tions over drj. This allows the spectra of photons and 
leptons to be presented in the form 



dN 

He 



H(r),E)dr). 



(69) 
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H( V , E y ) , o.u. 



1 r 




Figure 13: Function H{r], _E 7 ) at fixed energy of 7-rays, 
_E 7 = 1/2E+, calculated for a power-law distribution of pro- 
tons fp(Ep) oc E~ a with a = 2 and 2.5. 



Here 

where E is the energy of 7-rays or leptons, and $ is the 
energy distribution of the given type of particle. 

For photomeson interactions, it is useful to introduce 
the following characteristic energy of proton 

^=™ P c 2 ^ % )« 3.0.10 20 eV. (71) 

At energy E p = E*, the proton and a photon of energy 
kT can produce a pion through a head-on collision. 

The function H (77, E) at fixed energy E describes dis- 
tribution over 77. For a power-law distribution of protons, 
fp(Ep) oc E~ a , the function H(rj,E) has a maximum at 
v/Vo ~ 3 E/E*; the position of the maximum slightly de- 
pends on the power-law index a of the proton distribu- 
tion. The function H(rj,E) for 7-rays is shown in Fig. ^ 
at E-y — 0.5 E+ and two power-law indices, a — 2 and 2.5. 
At low energies the function H(rj, E) drops very quickly. 
A cutoff in the spectrum of protons results to faster de- 
crease of H (77, E) and to a shift of the point of maximum 
towards smaller 77. For electrons and neutrinos -ff (77, E) 
has similar behavior - a maximum at 77 <C 100 779 and 
strong decrease with increase of rj. 

To use the analytical parameterizations for distribu- 
tions of photons and electrons given by Eq.(27) and 
Eq.(31) and related Tables I, II, III, which are applicable 
for 77/770 < 100, the following condition should be ful- 
filled: 3E/E* <C 100. In this case the main contribution 
to the integral ( |69| ) comes from the region 77 <C 100 770, i.e. 
from events close to threshold. Therefore the obtained 
approximate analytical presentations allow calculations 
of distribution of particles in the energy range E < E*. 

Finally, let's discuss the production of photons and lep- 
tons at interactions of photons with a realistic distribu- 
tion of protons, namely a power-law with an exponential 



cutoff: 

f p (E p ) = AE; 2 exp (-;!£) • (72) 

where the normalization coefficient is determined from 
the condition 

00 

J E p J p (E p )dE p = 1 (73) 

lGcV 

In Figs. |l4] - |l7] we show the spectra of photons, elec- 
trons, and neutrinos produced in photomeson interac- 
tions calculated for 4 values of the cutoff energy in the 
proton spectrum E cut — 0.1 • E*, E*, 10 • E„ and 10 3 • E*, 
respectively. The case E cut = 10 3 • E* is almost identical 
to a pure power-law spectrum of protons. 

In Figs. |TJ| and we compare the spectra of electrons 
(and positrons) produced through the pair production 
channel with the spectra of electrons from the decay of 
photo- produced charged pions. 

Finally, in Fig. ^0], we show the steady-state spectra 
of cooled electrons. We assume that electrons are cooled 
via synchrotron radiation in the intergalactic magnetic 
field B — 1 /iG and inverse Compton (IC) scattering on 
the 2.7 K CMBR. Since the production spectrum of pair- 
produced electrons below E ~ 10 15 eV drops sharply (see 
Figs, [l^ and |l9|), the synchrotron and IC cooling (in the 
Thomson regime) leads to the formation of a standard 
E~ 2 type spectrum. This is clearly seen in Fig. 

In Fig. |2l] we show the spectra of synchrotron and 
IC radiation of secondary electrons produced via pair- 
production and photomeson production channels for a 
fixed magnetic field of B = 1 fiG, the temperature of 
CMBR T = 2.7 K and for 4 different cutoff energies in the 
proton spectrum E cut . Fig. |2l] a corresponds to the cutoff 
energy of protons E = 0.1E*. In this case the electrons 
are contributed mainly from pair-production process 
with a maximum in the energy distribution (E 2 dN/dE) 
at energy E ~ 10 15 eV. While synchrotron radiation of 
these electrons peaks at E 1 oc BE 2 ~ 10 6 eV, the max- 
imum of the IC radiation appears, because of the Klein- 
Nishina effect, at E 1 ~ E± ~ 10 14 - 10 15 eV. Note that 
although energy density of the magnetic field correspond- 
ing to B = 1/iG is B 2 /8tt « 4x 10~ 14 erg/cm 3 , i.e. an or- 
der of magnitude smaller than the energy density of 2.7 K 
CMBR, emissivity of the synchrotron and IC components 
are comparable. This is also a direct consequence of the 
reduction of the cross-section of IC scattering of 10 15 eV 
electrons in the Klein-Nishina regime. The second com- 
ponent of synchrotron radiation related to the electrons 
from photomeson processes peaks at much higher ener- 
gies, Ery no 10 12 eV, however its contribution is not sig- 
nificant because of suppression of the protons at energies 
above the threshold of photomeson reactions. The in- 
crease of the cutoff energy in the proton spectrum, E cut , 
leads to dramatic, orders of magnitude, increase of the 
emissivity of the synchrotron radiation of photomeson 
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electrons (see Figs. |2l|b,c,d). At the same time, because 
of the Klein-Nishina cross-section, only pair-produced 
electrons contribute to the IC radiation. Therefore the 
cutoff energy E cut does not have any impact on the IC 
spectrum and emissivity, as it is seen in Figs. ^l|. 

VI. SUMMARY 

We present simple analytical parametrizations for en- 
ergy distributions of photons, electrons, and neutrinos 
produced in interactions of relativistic protons with an 
isotropic monochromatic radiation field. The results on 
photomeson processes are obtained using numerical sim- 
ulations of proton-photon interactions based on the pub- 
lic available Monte-Carlo code SOPHIA. We also devel- 
oped a simple formalism for calculations of energy spec- 
tra of electrons and positron from the pair production 
(Bethe-Heitler) process based on the well-known differ- 
ential cross-section in the rest frame of the proton. The 
energy loss-rate of protons due to photomeson and pair- 
production processes in the 2.7 K CMBR calculated by 
integrating the energy distributions of the stable prod- 
ucts of interactions is in excellent agreement with re- 



sults of previous works based on direct calculations of 
energy losses (without intermediate stage of energy dis- 
tributions of secondaries). The analytical presentations 
of energy distributions of photons and leptons obtained 
in this paper provide a simple but accurate approach 
for calculations of broad-band energy spectra of gamma- 
rays, electrons and neutrinos in different astrophysical 
environments. 
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Figure 14: The production spectra (EdN/dE) of photons and electrons (left panel) and neutrinos (right panel) produced with 
energy distribution described by Eq.(Q) through the photomeson channel. The cutoff energy in the proton spectrum is assumed 
E cut = 0.LE*. 
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Figure 17: The same as in Fig. but for E cut = 1000 ■ £». 
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Figure 18: The production energy spectra of electrons and positrons produced through the channel of pair production (curve 
1) and positrons and electrons produced through the photomeson interactions of protons (curves 2 and 3, respectively). The 
proton spectrum is assumed in the form given by Eq.([72|) with cutoff energy at E cut = 0.1i5* (left panel) and iJ C ut = E t (right 
panel) . Note that the contribution of electrons (curve 3) in left panel appears below the low bound of y axis 
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Figure 19: The same as in Fig. ^ but for cutoff energies E cut = WE* (left panel) and E cut = 1000i5, (right panel) 
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Figure 20: The cooled spectra of electrons and positrons, N = iV+ + JV_. Dashed lines - electrons produced through the 
pair-production channel, dot-dashed lines - electrons produced through photomeson interactions. The sum of two contributions 
is shown by solid curves. The proton spectrum is given in the form of Eq.(|7^), with E cut — 0.1 ■ E* (a), E* (b), 10 • E* (c) and 
1000 -E, (d). 
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Figure 21: The synchrotron and IC spectra of cooled electrons. Dashed line corresponds to synchrotron (S) and IC radiation of 
electrons and positrons produced in the pair-production process, the dotted line corresponds to radiation of positrons (electrons) 
produced through photomeson interactions. Curve 1 (solid line) is the sum of these contributions. Curve 2 represents the 
spectrum of gamma-rays produced at decay of photoproduced 7r°-mesons. The proton spectrum is given in the form of Eq,(fr2|), 
with Ecut = O.IE* (a), E cu t = E, (b), E cut = 10 E* (c) and _B C ut = 10 E, (d). Magnetic field B = 1 fiG, temperature T = 2.7 K. 



